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ABSTRACT 

Although N-body studies of dark matter halos show that the density profiles, 
p(r), are not simple power-laws, the quantity p/cr^, where a{r) is the velocity 
dispersion, is in fact a featureless power-law over ~ 3 decades in radius. In 
the first part of the paper we demonstrate, using the semi-analytic Extended 
Secondary Infall Model (ESIM), that the nearly scale-free nature of p/o"^ is a 
robust feature of virialized halos in equilibrium. By examining the processes in 
common between numerical N-body and semi-analytic approaches, we argue that 
the scale-free nature of p/a^ cannot be the result of hierarchical merging, rather 
it must be an outcome of violent relaxation. The empirical results of the first 
part of the paper motivate the analytical work of the second part of the paper, 
where we use p/a^ oc r~" as an additional constraint in the isotropic Jeans 
equation of hydrostatic equilibrium. Our analysis shows that the constrained 
Jeans equation has different types of solutions, and in particular, it admits a 
unique "periodic" solution with a = 1.9444. We derive the analytic expression 
for this density profile, which asymptotes to inner and outer profiles of p ~ ^-o.is^ 
and p ~ r"^'^*^, respectively. 
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1. Introduction 



There is a broad consensus that gravity-driven cosmological evolution of coUisionless 
cold dark matter (CDM), starting from some realistic matter power spectrum, results in 
virialized objects whose spherically- averaged density profile is well represented by the NFW 
prescription (Navarro et al. 1997), or its close variants (Moore et al. 1998; Navarro et al. 

2004) . The main characteristics of the density profile are the relatively shallow density 
slope in the center, d\ogp/d\ogr ~ —1, and much steeper slope around the virial radius, 

—3. These general properties are well established, but the details are not. The value of 
the central slope, and the nature of the transition between inner and outer slopes is still a 
matter of debate. Currently, some advocate the original NFW formulation with constant 
asymptotic density slopes and a scale-radius (Diemand et al. 2005). Others argue for a more 
continuous change of the slope as a function of radius (Navarro et al. 2004; Merritt et al. 

2005) . Regardless of which density profile expression proves to best describe N-body halos, 
the underlying physics that drives the halos to have this shape is not yet fully understood. 

Several attempts have been made in the literature to explain the density profile shape. 
For example, Syer & White (1998) suggested that the NFW profile results from a series of 
repeated mergers; as the undigested cores of merging satellites accumulate at the center, 
the central r^^ cusp develops. However, Huss et al. (1999) performed simulations of isolated 
halos with varying degrees of thermal motion of dark matter particles. The final halos closely 
resembled the universal NFW shape, suggesting that hierarchical merging is not critical to 
the outcome. Moore et al. (1999) simulated halo formation using a matter power spectrum 
with suppressed small scale power. The final density profiles were still well described by the 
NFW fit. Some analytical studies suggest that secondary infall and accretion can determine 
the density profile shape (Avila-Rees et al. 1998). Barnes et al. (2005a) argue that the 
density profile, and in particular its radially changing slope, is the result of a mild aspect of 
the radial orbit instability. 

A different approach was taken by Taylor & Navarro (2001) (hereafter TNOl), and 
Hansen (2004). Instead of considering the radial run of the space density of N-body halos, 
they measured the quantity p(r) / a'^ {r) , which has the dimensionality of phase-space den- 
sity. In spherically-symmetric equilibrium halos, p/a^ is proportional to the coarse-grained 
phase-space density, a quantity distinct from the (fine-grained) phase-space density whose 
conservation is ensured by the coUisionless Boltzmann equation {e.g. Dehnen 2005). As our 
analysis does not hinge on the strict interpretation of the quantity p/cr^, we will continue to 
call it phase-space density, for brevity. 

TNOl found that p/a^ of their three halos (circular velocities 160-200 kms~^) is very 
well approximated by a single power-law, from 10~^"^ virial radii to just beyond the virial 
radius. This finding is unexpected because the density profile, p{r) undergoes a considerable 
slope change in the same radial interval. The scale-free nature of the phase-space density 
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implies that the double logarithmic slope of the velocity dispersion changes in such a way so 
as to offset the change in the density profile slope. The power-law index a, where p/a^oa r~", 
was found to be very similar for all the three halos studied by TNOl; a — 1.875. 

Does the power-law nature of p/a^ characterize all equilibrium N-body halos, or are the 
halos studied by TNOl somehow special? One can get a rough answer to this question without 
analyzing many individual simulation halos. Given the form of a "universal" empirical fit to 
N-body halos and assuming that the halos are in equilibrium, one can calculate the phase- 
space density as a function of radius using the Jeans equation. Barnes et al. (2005b) use two 
types of empirical fits to N-body profiles (Navarro et al. 1997, 2004) to demonstrate that 
both produce nearly scale-free phase-space density profiles for ~ 3 decades in radius, with a 
slope close to 1.9. They conclude that most N-body halos should have scale-free phase-space 
density profiles. This conclusion is further supported by a recent study of equilibrium halos 
(Kazantzidis et al. 2004). The halo that most resembles results of N-body simulations has 
scale- free p/a^ and index a 1.95. 

The scale-free nature of p/a^ represents a novel way of looking at the properties of 
halos, and as such, should be explored further. If this property is "universal" , it amounts 
to a hitherto unrecognized constraint on the shape of the density profiles. In this paper we 
address questions related to the p/a^ profiles of coUisionless halos. Is the scale-free nature of 
p/a^ (i) a result of initial conditions, (ii) an outcome of evolution and virialization, or (ni) a 
consequence of the final equilibrium state of halos? In a related work, Barnes et al. (2005b) 
address the third question. They demonstrate that the condition of hydrostatic equilibrium^ 
alone, as quantified by the Jeans equation, does not require p/a^ to be a power-law in radius. 

The first two questions can be addressed by studying halo formation techniques different 
from N-body simulations. Semi-analytical methods based on the secondary infall paradigm 
are an example of such techniques. Halos forming in N-body simulations are subject to a 
complex set of dynamical processes. Particles can exchange energy and momentum through 
violent relaxation and two-body encounters (even though the latter are suppressed in cos- 
mological simulations). Additionally, phase- mixing can decrease phase-space density, while 
preserving particles' energies. These processes are brought about by halo collapse, merging, 
dynamical friction, etc. The analytical methods do not incorporate such a wide range of 
processes. The method we use here is based on Ryden & Gunn (1987), and was extended 
by Williams et al. (2004b); we call it ESIM, the Extended Secondary Infall Model. As we 
describe later in more detail, ESIM halos are spherically symmetric with one spatial and 
two velocity (radial and tangential) phase-space dimensions. They virialize through violent 
relaxation, which changes the shells' energies, but conserves the angular and radial momenta 
of each shell throughout the collapse. There are no mergers. 



^We use the phrase 'hydrostatic equihbrium' because coUisionless dark matter can be thought of as fluid, 
with pcr^ acting as pressure (Shapiro et al. 2004). 
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Given such different formation scenarios, it is not surprising that the final properties — 
density distributions, etc. — of N-body and ESIM halos tend to be different. (See Barnes 
et al. (2005a) for a discussion of why the two types of density profiles differ.) However, 
if specific final properties are present in N-body and ESIM halos, that strongly suggests 
that the property is the result of shared physics. This is our premise in addressing question 
(ii). Specifically, wc find that the phase-space density of our halos is well represented by a 
power-law (Section 2), a result similar to the one reached through the study of numerically 
generated halos. Question (i) is also addressed in Section 2, where we vary input conditions 
of ESIM halos. Again, the power-law nature of p/a^ does not seem to be affected by these 
indicating that the answer to question (i) is negative. 

In Sections 3 we address the behavior of index a. One of the main properties of the viri- 
alized halos is that they are in equilibrium, so the Jeans equation is important in determining 
the structure of halos. Therefore it makes sense to ask if the Jeans equation has a preferred 
value for a, or, if all values are equally likely from the point of view of the Jeans equation. 
Our analysis of the second derivative of the Jeans equation and its various solutions leads to 
some interesting insights (also see Appendix A). Appendix B derives the value of a for halos 
formed in a highly idealized secondary infall model (SIM), and compares it to a's obtained 
in ESIM halos of Section 2. 

2. Extended Secondary Infall Model (ESIM) Halos 

2.1. Summciry of the Method 

The collisionless dark matter halos in the present study arc formed using an analytical 
method originally introduced by Ryden & Gunn (1987). The method, including minor 
extensions, is fully described in Williams et al. (2004b). In summary, halos are spherically 
symmetric and isolated. We use a fiat, matter dominated geometry with h — 0.5. The initial 
power spectra is a standard CDM Harrison-Zeldovich power spectrum normalized such that 
Q°-^(T8 = 0.5. Evolution of the power spectrum is dictated by the transfer function given 
by Bardeen et al. (1986) in their Appendix G. The density profile of a proto-halo is taken 
to be the profile of a peak in a density field described by this power spectrum, P{k) as is 
illustrated in WiUiams et al. (2004b), Figure 9. 

In addition to the primary peak, proto-halos have superimposed secondary perturba- 
tions, whose properties are derived from the same P{k). In a real halo these secondary 
perturbations generate random motions of dark matter particles, both radial and tangential. 
In ESIM halos, the RMS amplitude of random velocities (as a function of radius) are calcu- 
lated from P{k), as described in Ryden & Gunn (1987), and in the Appendix of Williams 
et al. (2004b). For each shell the magnitude of the velocity perturbation is chosen from a 
Maxwell-Boltzmann distribution with the above RMS; the direction is chosen randomly. It 
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has been recently pointed out that the Maxwell-Boltzmann distribution does not describe 
the particle velocities of equilibrium halos (Kazantzidis et al. 2004; Hansen et al. 2005). 
However, the Maxwell-Boltzmann distribution is sufficient for our purposes, because we use 
it to describe the initial pre-coUapse, not the final equilibrium halos. Even if the Maxwell- 
Boltzmann assumption is not strictly correct for the proto halos, it suffices for our purposes, 
because we use ESIM halos as a simplified tool only; we do not aim to reproduce real halos. 

When a particular shell reaches its maximum turn-around radius the perturbation in 
radial and tangential velocity is added to the existing radial infall velocity. In a real halo 
these additional random velocities result in dark matter particles having elliptical, precessing 
orbits. In ESIM halos, the random velocities transform the initially thin shells into thick, 
overlapping shells, each with its own apo- and peri-center. A key assumption in ESIM 
calculations is that the collapse proceeds adiabatically, with each shell conserving its angular 
and radial momentum. Note that each shell's secondary velocity perturbation is oriented 
randomly, so that the halo as a whole has zero net angular momentum, and no net rotation. 
As the collapse proceeds, the mass distribution in the halo changes, and the potential and 
shells' energies are recalculated after each step in the collapse. This is violent relaxation. To 
sum up, ESIM halos are formed via smooth accretion of somewhat lumpy material, with no 
major mergers. 

Our method of halo formation has considerable flexibility with direct control over the 
parameter space of initial conditions. For example, we can change the shape of the primary 
density peak by smoothing or altering the matter power spectrum. Both the shape of the 
primary peak and the amplitude of the secondary perturbations are based on P{k), but 
we can alter the shape of the primary peak and the amplitude of secondary perturbations 
independently, which gives us another degree of freedom in exploring the resulting final halo 
profiles. For example, we can change the secondary perturbations — the overall amplitude and 
radial dependence — while keeping the primary peak unaltered. Even though these changes 
need not correspond to any astrophysically plausible scenarios, the evolution of ESIM halos 
proceeds in a physically consistent way. 

2.2. ESIM Halo Calculation Results 

We generate three sets of halos, with 'galaxy', 'group' and 'cluster' scale masses. To get 
different ranges of halo mass we filter the mass power spectrum, prior to normalization, on 
comoving scales rh of 0.05, 0.35 and 0.7 Mpc, respectively. Larger Vh mean that the collapse 
will start later, and more massive halos will form as a result. Within each mass set we study 
the effect of the secondary velocity perturbations on the final halo profiles. To that end, we 
form halos where the magnitudes of the secondary velocity perturbations are scaled by a pa- 
rameter, hereafter called the perturbation amplitude scaling parameter, or scaling parameter 
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for short. In this Section the scahng parameter is altered uniformly for a given halo, but in 
Section 2.3 we radially vary it within the halos. A scaling parameter of zero corresponds to 
no secondary perturbations, in which case the velocities are purely radial. A scaling param- 
eter of 1 corresponds to the original amplitude of secondary perturbations derived from the 
same power spectrum that determined the shape of the primary density peak, while setting 
the scaling parameter to 2 doubles the amplitude of all secondary perturbations. 

The perturbation amplitude scaling parameter is varied from 0.4 to 2.5 in increments of 
0.05. Because each shell is assigned a velocity perturbation from a distribution extending to 
is a chance that the shell's kinetic energy will exceed its potential energy. 
This is especially likely if the scaling parameter is large, and the shell happens to be far 
from center. In this case, we assume that the shell is instantaneously lost from the halo 
and is no longer considered in the potential calculations or the kinematic analysis. Each 
shell's potential energy is calculated based on the mass distribution due to a finite number 
of shells, with the outermost shells extending well beyond the viriahzed portion of the halo 
at the present epoch. The total number of shells, and hence the zero-point of the potential, 
is fixed such that halos computed using scaling parameters less than 1.5 do not lose any mass 
during the collapse, while halos with a scaling parameter of 2 lose at most 25% of their final 
mass. 

Each ESIM halo is unique since the secondary velocities for each shell are picked ran- 
domly from a distribution. To reduce the dispersion in the final halo properties the halos 
presented here are averages of ten individual halos. To verify the reliability of the results two 
sets of averaged halos were computed for each scaling parameter. In all cases our analysis 
applies to the portion of each halo between 10^^'^ and 1 virial radius, r2oo; the radius at which 
the average enclosed halo density is 200 times the critical, Po = 2.77 x IO^^H'^Mq Mpc^^. As 
mentioned earlier, we use a standard cold dark matter model. While this model is not the 
currently favored one, the particular cosmological model used has a minimal effect on the 
final halo properties of interest. 

In the following sections we describe the various properties of the final virialized ESIM 
halos. In Sections 2.2.1-2.2.4 we discuss in detail the properties of galaxy mass halos. Since 
the qualitative properties and the general trends are similar for galaxy, group, and cluster 
mass halos, we do not discuss the latter two specifically, but present a summary of their 
properties in Section 2.2.5. 

2.2.1. Density profiles 

Figure 1 shows four halo density profiles, expressed in units of critical density Po, and 
multiplied by {r /r2QoV^''* i where — 7/it is slope of the best power-law fit to the density profile 
between 10~^"^r/r2oo and r/r2oo- This is done in order to accentuate features in the density 
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profiles. All halos were generated using the galaxy-mass P{k) filtering scale rh — O.OSMpc. 
The standard case, a halo with a scahng parameter of 1, is shown in panel C. Panels A, B 
and D show halos with scahng parameter of 0.4, 0.6 and 2, respectively. 

As in Williams et al. (2004b), the ESIM halos differ from those in N-body simulations. 
For convenience, NFW density profiles, p{r) oc (r/rs)(l -\-r/rsY (Navarro et al. 1997), with 
concentration parameter (c = r2oo/'"s) appropriate for that halo mass, and arbitrary normal- 
ization are over-plotted as long dash fines. ESIM density profiles are better approximated 
by a single power-law between r = 10~^'^r2oo and r2oo than are NFW profiles. 

The effects of changing the amplitude of secondary perturbations on the mass distri- 
bution in the halos can be seen in the four panels of Figure 1, and is summarized in panel 
A of Figure 2, which plots 7/jt, the fitted density profile slopes, vs. scaling parameter. As 
expected, the halos with larger secondary perturbations, and hence more angular momentum 
have shallower density profiles. (Note that the steepest central cusp slope of ESIM halos is 
7 = 2, as explained in Section 5.1 of Williams et al. (2004b).) To further characterize the 
density profiles we indicate in Figure 1 where the slope, = — ((ilnp)/((ilnr), attains the 
values of 7=1, 2, and 3, with dotted, short-dash and long dash fines. We denote these radn 
by r^=i, r^=2-i and r^=3. The dependence of r.^=i^2,3 with respect to the scaling parameter 
is shown in panel D of Figure 2. The values of r-y=i,2,3 increase with the scaling parameter, 
consistent with density profiles becoming shallower as the angular momentum content of the 
halos increase. 



2.2.2. Velocity dispersion and velocity anisotropy 

Since many shells overlap at any given radius, the total velocity dispersion profile is 
calculated as the mass- weighted and time-weighted RMS value of the radial and tangential 
velocity dispersions of all the shells that exist at that radius. Figure 3 shows the velocity 

dispersion profiles for the halos presented in Figure 1. The values of r^=i, r^=2, and r-y=3 are 
marked as dotted, short-dash and long-dash lines. Figure 3 also shows the velocity dispersion 
anisotropy, (3, defined as usual, through the radial and tangential velocity dispersions (Binney 
& Tremaine 1987), 



(3 = 1- 




Except for halos generated using small scaling parameter values, < 1, the anisotropy param- 
eter stays close to between ~ 10~^ and ~ 1 virial radii. For the standard halo, the average 
velocity anisotropy parameter is P ~ 0.1. Beyond the virial radius, the velocity ellipsoid 
parameter increases to 0.5 — 1 as the orbits become more radial. In general, a halo's average 
(3 increases with decreasing scaling parameter: smaller secondary perturbations imply less 
angular momentum, and hence more radial orbits. In the limit of perturbation amplitude 
scaling parameter of 0, purely radial orbits will result. 
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In Figure 4, we show the dependence of the velocity dispersion anisotropy with respect 
to the density slope, 7/it. This relation has been presented for N-Body simulations by 
Hansen & Moore (2004). They have found that there is a tight correlation between /3 
and 7 irrespective of the considered models. ESIM halos do not display this correlation. 
One possible explanation for this difference between ESIM and N-Body halos is that ESIM 
formalism docs not allow the onset of the radial orbit instability, which may be key in forming 
/? — 7 correlation in N-Body halos (Barnes et al. 2005a). 



2.2.3. Phase-space density 

The phase-space density profiles of the halos displayed in Figm^c 1 are shown in Figure 5. 
ESIM halos appear to have nearly scale- free distribution of p/a^ extending over 2 — 4 decades 
in radius, interior to the virial radius. Aside from the very central and outer regions (which 
we will return to shortly) p/a^ arc fit by power-laws better than the corresponding density 
profiles. For example, between 10^^'^ and 10~^'^ virial radii the density profile of the halo 
with the scaling parameter of 1 (panel [C]) is becoming shallow rapidly, but p/cr^ is still well 
fit by the same slope that characterizes larger radii. 

To understand why p/a^ deviates from a power-law at large and small radii we have to 
consider the equilibrium state of halos. Binney & Tremaine (1987), under the assumptions 
of spherical symmetry and a non-rotating system, rewrite the Jeans equation as 

l d{pal) ^ ^pal _ dcj) 
p dr r dr 

In Figure 6, the LHS of equation 2, i.e. the pressure gradient plus the anisotropy term, 
is plotted at each radius as points, while the RHS, the potential gradient, is shown as the 
smooth curve. Hydrostatic equilibrium is satisfied when the points follow the curve. This 
is the case in panels [A], [B], and [C]; in panel [D] there is a lot of scatter in the points 
that represent the LHS of eq. 2, and so the equilibrium condition is not well satisfied. This 
apparently results in p/a^ being only marginally well approximated by a power- law. In 
panels [A]-[C], the hydrostatic equilibrium breaks down outside of the virial radius, which 
apparently leads to a break in the scale-free behavior of p/a^, at the same radius. 

Figure 6 also helps to explain why p/a^ breaks away from a power-law at small radii. 
At small radii, the potential gradient scales as This is due to the presence of a compact 
inert core at the center of every ESIM halo. As described in Williams et al. (2004b) (Sections 
4.1, A.l, A. 5), ESIM halos have finite mass resolution. At the center of each halo the smallest 
resolvable mass is represented by an inert core, which provides the initial 'seed' for ESIM 
calculations. The inert core does not evolve, i.e. its mass distribution does not change as 
the halo collapses. Because the part of the halo where the core dominates is not subject to 
the same dynamical evolution as the rest of the halo, that part should not be considered in 



8 



the analysis. In panels [A]-[D] of Figure 6 the regime of core domination begins roughly at 
10~^, 10~^'^, 10~^, and 10~^"^ virial radii, respectively. The corresponding panels of Figure 5 
show that these radii coincide with the break in the scale- free nature of pjo^ . 

Excluding the innermost and the outermost regions of halos, the remaining portions 
are well represented by power-laws (this is further illustrated in Section 2.2.5 and Figure 9 
below). Figure 5 demonstrates that the scale-free nature of pja^ is resihent: it persists 
despite the changes in the magnitude of the secondary perturbations, at least in the range 
0.5-1.5. However, the power-law index, a = —{d\iap/a^)/{dlnr), is not constant; in general, 
a fit increases with increasing scaling parameter (see panel B of Figure 2). 

2.2.4- Halo mass 

Panel D of Figure 2 shows that the mass of virialized halos, for a given power spectrum 
filtering length (r^ = 0.05Mpc for these halos), decreases steadily as the perturbation ampli- 
tude scaling parameter increases. Large secondary perturbation amplitude implies a more 
extended mass distribution, and smaller virial masses. Additionally, if the shells' energy 
reaches the threshold to cause mass loss, the final mass is reduced further. However, this 
effect is negligible for most halos: there is no mass loss during formation for halos with the 
scahng parameter less than about 1.5. 

2.2.5. ESIM halos of galaxy groups and cluster masses. 

So far, we have discussed how the properties of galaxy-mass halos {vh = 0.05Mpc) 
change with the changing amplitude of secondary perturbations. Group and cluster-mass 
halos, generated using power spectrum filtering lengths of r/j = 0.35 and 0.7Mpc, respectively, 
behave similarly. The overall trends for these two cases are summarized in Figures 7 and 8, 
which are similar to Figure 2. 

Two examples arc shown in Figure 9, for Vh = 0.35Mpc (panels A through E) and 
r/i = 0.7Mpc (panels F through J); both halos have scaling parameter of 1. It is apparent 
that larger r/j values result in density profiles that are not as well described by a single 
power-law, as smaller rh halos. This is further illustrated in panels C of Figure 7 and 8; 
'"7=1 ) '"7=2 

and r-y=3 are closer together than in galaxy-mass halos, implying that the density 
slope changes more rapidly with radius. 

Despite the non power- law behavior of density, the phase-space density, p/cr^, remains a 
power- law for at least 2-3 decades in radius interior to r2oo (panels D and I of Figure 9). These 
panels also show that the scale-free nature of p/a^ breaks down interior to 10~^-^ — 10~^, 
where the inert core begins to dominate the potential (panels E and J) . The same behaviour 
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was noted in the case of galaxy-mass halos, in Section 2.2.3. The best- fit slope afu increases 
with the scaling parameter for a fixed Vh (panels B of Figure 2-8), but gets consistently 
smaller with increasing r/j. 

Figure 10 depicts the dependence of the p/a^ power-law slope, afu on the virial mass 
for all three sets of halos: Vh — 0.05 Mpc (triangles), — 0.35 Mpc (squares), and Vh — 0.7 
Mpc (pentagons). Filled symbols show afu obtained by fitting a power-law to the data 
within our standard radial range, from 10~^"^ to 1 virial radii. To test the dependence of 
the specific choice of the radial range we also use a reduced range of 10~^'^ to lO""'^ virial 
radii (empty symbols). The afu values are similar in both cases. The main conclusion from 
Figure 10 is that aju is not uniquely determined by any single halo property, such as mass, 
perturbation amplitude scaling parameter, or Vh- 

2.3. Testing the Resiliency of the p/a^ Power- law 

Since the halos considered in previous sections exhibit a near-perfect power-law in p/a^, 
we wish to determine what would be required to break this resilient behavior. To that end, we 
generate halos with different scaling parameters assigned to blocks of shells. These changes 
are not meant to mimic realistic situations, but rather to push our simulations to unphysical 
extremes. In several trials we increased or decreased the scaling parameter linearly with 
radius; in other trials, we assigned the inner 5-15% of radii a larger perturbation amplitude 
scaling parameter of 5, 20, or 50, with smaller scaling parameters assigned to outer shells. 

The halos with the scaling parameter assigned linearly in radius (with no large disconti- 
nuities) display the power-law in p/a^. However, halos in the trials with the 'excited sections' 
(blocks of shells with large scahng parameters) had a, p/a^ profile that showed deviations 
from a power-law behavior. Figure 11 shows the specific halo where all the shells were as- 
signed a scaling parameter of 0.9, except for those with an initial comoving radius between 
0.59 and 0.79 Mpc, which were assigned a scaling parameter of 5. The last shell to collapse in 
this evolution had a comoving radius of 1.92 Mpc. Visual inspection of panels A and B show 
that, even though these halos have been subjected to radially discontinuous changes in the 
amplitude of the secondary perturbations, the final density and velocity dispersion profiles 
still appear well behaved, but the deviations of p/a^ from a power-law are evident. This and 
other similar trials illustrate that one must adopt extreme measures — large discontinuities 
in the scaling parameter — in order to break the scale-free nature of the phase-space density 
profile. 
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2.4. Conclusions from ESIM Halos 



The ESIM formalism provides us with a simphfied, but physically self- consistent way 
of generating coUisionless dark matter halos. It is very different from N-body simulations. 
In the preceding sections we described ESIM halo properties in detail, pointing out that 
the density profiles are generally different from those produced in N-body simulations. By 
changing the initial conditions we produced a range of density profile types: by filtering 
the mass power spectrum on small scales {rh — O.OSMpc) to simulate galaxy-mass halo 
formation, we ended up with density profiles that showed virtually no change in slope over 
2-3 decades in radius inside of the virial radius (Figure 1). Larger filtering scales {rh — 0.35 
and 0.7 Mpc) aimed to produce group and cluster-mass halos resulted in more curved density 
profiles (Figure 9). 

The most important conclusion of Section 2 concerns the phase-space density profile, 
p/(T^. Despite the differences in methods and density profile shapes, one feature was shared 
by all ESIM halos and N-body simulated halos: the power- law nature of the p/a^ profile. We 
have shown that this feature of ESIM halos is very resilient, and one must concoct rather 
extreme cases in order to end up with p/cx^ that shows deviations from a single power- 
law (Section 2.3). There are no models in astrophysics literature that provide a physical 
explanation for this feature of virialized halos. We argue that the scale-free nature p/a^ 
cannot be uniquely characteristic of hierarchical merging because ESIM halos do not undergo 
mergers. We speculate that it is a generic consequence of the one dynamical process which 
is common to both ESIM and N-body halo formation — violent relaxation, or, global changes 
in the potential, and the accompanying changes in particle energies brought about by the 
collapse. 

Our results concerning a, the double logarithmic slope of p/a^ vs. radius, are different 
from those of TNOl. We find a range of values from 1.6 to 2.2, while TNOl find that the same 
a = 1.875 described all three of their N-body halos. It is possible that the latter finding 
is due to a small number of halos studied. In our study, Ofu was not uniquely determined 
by any single property of halos, such as mass, epoch of formation, or amplitude of random 
velocities. 



3. Jeans equation analysis 

3.1. The second derivative of the Jeans equation 

The findings of the preceding section, summarized in Section 2.4, motivate our next step. 
ESIM halos are spherically symmetric, in equilibrium, have zero net angular momentum, and 
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have approximately isotropic velocity dispersion ellipsoids, P ^ 0.^ These halo properties 
lead directly to the appropriate version of the Jeans equation, and its first derivative: 



GM{< r) 



(3) 



dr 



which states that hales are supported against collapse by "pressure," provided by the velocity 
dispersion of dark matter particles. Because halos are assumed to have isotropic velocity 
dispersions, a — ar — ere — Ufj,. In Section 2 we have established empirically that ESIM 
halos exhibit scale-free p/ behavior. In this paper we will not address the question of why 
phase-space density radial profiles are power-laws. Instead, we will take this property as an 
additional constraint: 



where ^O; ctq; and a are constants for any given halo. Eliminating cr, and using the 
dimensionless variables, x — r/ro and y = p/po, we can now write down the derivative of 
equation 3, 



This is also equation 4 of TNOl. The constant n = 47rGpo'^o/'^o specifies the normalization 
of halos. Exponent a is constant for any given halo, but we leave its value unspecified. 

In the rest of this Section we develop a framework for studying halos that can be 
described by equation 5. However, unlike the earlier Jeans equation work our main equation 
will not be equation 5 itself, but its first derivative, which is also the second derivative of 
the Jeans equation: 



where we have defined logarithmic derivatives of the density slope, 7 = —dlny /dlnx, as 
7' = dj/dlnx, and 7" = d'-f'/dlnx. Equation 6 has several advantages over equation 5. 
It docs not depend on the normalization constant k, and has no explicit dependence on 
X, in fact, equation 6 depends only on the exponent a, and the local density slope, 7(2;). 
Equation 6 looks like an algebraic equation — its partially factorized form will prove useful 
in exploring the various types of solutions of the Jeans equation. 



^Looking at Figure 3 zero anisotropy is arguably the least justified of our assumptions. N-body generated 
halos have velocity dispersions that are nearly isotropic at the center, and become radially anisotropic further 
out. Strictly speaking, the analysis in this Section applies only to the portions of halos that are isotropic. We 
will extend our analysis to arbitrary anisotropy in the future work. As a place to start, /3 w is a reasonable 
approximation. 

^This equation corrects a sign error of eq.(2.2) of Williams et al. (2004a). Note that equation 6, together 
with some interesting analysis also appears in Dehnen & McLaughlin (2005). 




(4) 




(5) 



(2q; + 7 - 6){2a - 57)(|[a - 7] + 1) = 87 '(80; - 57 - 5) + I57", ^ 



(6) 
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Our goal in the following is to explore solutions of the Jeans equation in general, paying 
particular attention to the role of a. Specifically, we want to know if the Jeans equation 
admits any special a values. 

3.2. The q; — 7 plane 

A halo is characterized by a relation between 7 and x, which can be obtained by assuming 
a value for a, specifying initial conditions for 7 and 7', and then integrating equation 6. An 
infinite variety of halos can be generated in this fashion, so a straightforward exploration of 
the parameter space is probably not the most insightful one. The form of equation 6, namely 
that it looks like a factorized algebraic equation of two variables, a and 7, suggests that it 
is instructive to consider the a — 7 plane, shown in Figure 12. Assuming a is constant for 
any given halo, a halo density profile can be represented by a vertical straight line segment. 

There are four factors in equation 6, three on the LHS, and one on the RHS. By sep- 
arately setting these to zero we obtain the four thick lines in the a — 7 plane. This plot 
provides a framework for understanding and classifying different types of solutions of the 
Jeans equation. 

Consider the line 2q; — 57 = (lower thick solid line). If part of a halo is described by 
a constant 7 = (2/5)q; for some range of radii, then the LHS of equation 3 can be expressed 
as d\\i{pa^) / dXnr (x 2a — 57, implying that the pressure gradient vanishes in that part of 
the halo. An example of this is the central region of the critical solution identified by TNOl, 
where 7 = 0.75, corresponding to a = 1.875. In general, halos of any a whose inner profile 
slope is described by 2q; — 57 = have vanishing pressure gradient in the center. 

The role of terms 2a + 7 — 6 and 8a — 57 — 5 is addressed in Sections 3.3 and 3.4 
respectively. The term |[a — 7] + 1 is briefly discussed in Section A, and all four terms of 
equation 6 are important in the discussion of Section 3.4. 

We note that halos with radially increasing density cannot exist; Tremaine et al. (1994) 
showed that if the distribution function depends solely on energy, as is the case with our 
isotropic halos, then 7 < is unphysical. 

3.3. Power-law density profiles 

One of the factors on the LHS of equation 6 is 2q; + 7 — 6; 2a + 7 — 6 = is shown as 
the long-dash line in Figure 12. 7pz = 6 — 2a represent power-law density proflle solutions 
to the Jeans equation (TNOl, Hansen (2004)), i.e. if 7' = 7" = then 7pz=const for all 
X. A power-law density profile is represented by a point in Figure 12. For example, the 
singular Isothermal Sphere is a viable solution of equation 6, characterized by a = 7 = 2 



13 



and constant a throughout the halo. 

Because equation 6 docs not distinguish between the three factors on the LHS, one 
might think that each of these three should produce power-law density profiles. This is not 

2 

so. Equation 5 can be rewritten as (2a — 57) (| [a — 7] + 1) — 57' = — a;^~3" y^l^ . If the 
slope 7 is given by 2a — 57 = 0, or | [a — 7] + 1 = 0, and 7 ' = because we are dealing with 
a power-law density profile, then the RHS becomes zero, which is unphysical. (Recall that 
y is scaled density.) Therefore, for a given a there is a unique power-law density solution, 
7pi = 6 - 2a. 

The one-to-one correspondence between 7^; and a can be tested using ESIM halos. 
Galaxy-mass {jh — 0.05 Mpc) ESIM halos presented in Section 2.2 have nearly power-law 
density profiles. Of these, halos with the perturbation amplitude scaling parameters of 1 or 
larger have nearly isotropic velocity dispersion ellipsoids. For these halos the relation between 
the fitted slopes of the density and phase- space density profiles should follow 2afit+jfit — 6 = 
0. In Figure 12 we plot afu and 'jfu for halos with scaling parameters between 1 and 2.5 
(crosses). There is good agreement between ESIM halo data and predictions. 

3.4. The unique periodic solution of Equation 6 

Equation 6 is a non-linear damped oscillator equation for 7, with the 7' term being 
analogous to a dissipation term. Equation 6 therefore suggests that the Jean's equation 
may admit solutions where the slope of the density profile oscillates with radius, becoming 
alternately steeper and shallower. If equation 6 is in a regime where these oscillations are 
damped, then the slope 7 may oscillate, with the amplitude of the oscillations decaying with 
increasing radii. Asymptotically, this behavior would lead to a density profile that has a 
power-law fall-off at large radii. Instead, if equation 6 is in a regime where the oscillations 
can amplify, then the slope may run away to large values, leading to a density profile that 
becomes increasingly steep at large radii. 

The coefficient of the 7 ' term in equation 6 controls the "dissipation" in the oscillations 
of the density profile slope 7. If this coefficient were always positive (negative), it would 
damp (amplify) the amplitude of 7 oscillations with Inx.^ If the coefficient is zero, the 
amplitude of oscillations would not evolve, resulting in perpetually periodic oscillations. In 
the case of equation 6, the damping term is non-linear, with its sign changing as a function 
of Inx. Does a periodic solution exist in such a situation? 



''We stress that the halos arc in cquihbrium, and the use of the word "oscinations" docs not mean that 
the halos are changing in time: the changes are in the density profile slope, as a function of logarithmic 
radius. 
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Let us use Figure 12 as a guide. Notice that the point of intersection of So; — 57 — 5 = 
and 2(3; + 7 — 6 = hues is exactly half-way between 2q; — 57 = and |[q; — 7] + 1 = 0. The 
corresponding a — 35/18 = 1.9444, and with that equation 6 can be rewritten as 



where B — '^pi — 'y, and -jpi — 2.111 for a — 1.9444. Equation 7 says that when B — 0, 
then 7" = 0, and hence 7' attains a local minimum or maximum. So, every time 7 = 2.111, 
7 ' attains an extremum. This statement can be true only if 7 varies periodically with In x. 
Getting back to the coefficient of the damping term: when a — 1.9444, the value of Set— 57— 5 
is, on average, zero. 

Thus, a — 1.9444 results in periodic solutions, and Figure 12 tells us that this a value 
is unique. The amplitude of the wave, i.e. the minimum and maximum values of 7, and its 

period of oscillations, i.e. the radial span in the halo between steepest and shallowest slopes, 
depend on the initial conditions. The maximum amplitude is attained for initial 7' = 0, 
and initial 7 somewhat larger than (2/5)a (lower thick solid line in Figure 12), or somewhat 
smaller than a + 1.5 (upper thick solid line in Figure 12). As the initial 7 is set closer to 
Ipi = G — 2a = 19/9 ~ 2.111 the amplitude of oscillations decreases. For the exact jpi we 
get a power-law density profile, as already noted in Section 3.3. The periodic behavior is 
illustrated in Figure 13 which shows the hmit cycle of equation 6, as long-dash and dotted 
line closed curves. For both of these a — 1.9444 and initial 7' = 0, but initial 7j„'s are 
different (0.7788 and 2.7111). 

The shape of the upper portion of the a — 1.9444, 7j„ Ki (2/5)a (=0.7788) curve in 
Figure 13 (long-dash line) suggests that the corresponding solution of the Jeans equation 
has a simple analytical form — a parabola. In fact, equation 6 does admit the following 
solution: 

7 =T„-5(7-7o)^ 70 = ^5(55 - 12a), 7„ = f+ ^ ^(6. + 1) ' 
Equations 8 allow us to express the density profile slope as a function of radius in the halo: 

7=(67j^/' tanh[(7:„/6)^/Mnx]+7o. (9) 

This solution is vahd exactly only for the halo with a — 1.9444 and asymptotic density 
slopes, 7j„ = 0.7778 and ^out — 3.444, given by the intersection of the a — 1.9444 hne with 
the upper and lower thick solid lines in Figure 12, respectively. For a's within ~ 1% of 
1.9444 equations 8-9 are a good approximation to the actual solutions of the Jeans equation. 
For other values, the approximation is less valid. For example, for a = 1.875 the above 
expression is shown as the thin long-short-dash line in Figure 14. For reference, we plot the 
value of 7o, the slope where the density slope changes most rapidly, as a thin solid line in 
Figure 12. As expected, for a = 1.9444, 70 = 7^; = 2.111. 
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Finally, we note that a — 1.9444 periodic solution value applies to the case of isotropic 
velocity dispersion, (5 — If a halo is characterized by a larger, but constant anisotropy, 
the periodic a value decreases. In the hmiting case of radial orbits and (3 — 1, a — 1.7222. 
Dark matter halos with non-zero anisotropy will be discussed in a future paper. 

3.5. Two classes of solutions 

It is clear from Figure 13 that a — 1.9444 is the watershed value that separates two 
different types of Jeans equation solutions: an example of a > 1.9444 is the short-dash hne 
spiral curve that eventually settles to a constant finite value of 7; an example of a < 1.9444 
is the solid line which shows a rapid increase of 7 to very large values. The damping term 
line, 8a — 57 — 5 = in Figure 12 explains this behavior. Solutions with a < 1.9444 tend 
to have a negative average damping term (and density profiles that steepen without limit), 
while solutions with a > 1.9444 tend to experience positive average damping (and result in 
density profiles that asymptote to power- laws, albeit at very large radii). We will call these 
the negative and positive damping solutions, respectively. 

Let us consider negative damping solutions and, for concreteness, set a = 1.875. We 
integrate equation 6 from the center of a halo outward, i.e. from small to large values 
of Inx. We start from given by 2q; — 67 = 0, and 7' slightly greater than 0, which 
ensures that the center of the halo will have constant density profile slope, 7j„ = (2/5) a, and 
hence vanishing pressure gradient (Section 3.2). As \nx increases, 7 increases without limit, 
because of the negative average damping term. The solid thick hne in Figure 14 shows the 
resulting profile. The same profile is also shown in Figure 13 with the same line type. (Note 
that the horizontal normalization is not important, and was adjusted for easier comparison 
between various results.) 

To illustrate the behavior of positive damping term solutions we set a — 1.975. As 

before, we choose 7j„ = (2/5)q;, and initial 7' slightly greater than 0. This solution is 
plotted as the short-dash line in Figures 13 and 14. It does not attain steep density slopes 
in the outer halo; 7 climbs to about 2.45, then oscillates, with decreasing amplitude, and 
eventually settles to the asymptotic value given by 2Q; + 7p; — 6 = 0. It is the positive average 
damping term in equation 6 that forces the oscillations to a constant asymptotic 7^;. The 
halo represented as the short-dash line has an outer density slope 7 < 3, implying that its 
mass does not converge. However, this is not a problem because in a real halo the profile 
pictured here would not apply at arbitrarily large radii where the assumptions of the Jeans 
equation will break down, including the condition of equilibrium itself. A real halo will have 
a steep outer slope, and convergent total mass. 

Thus, the behavior of 7 vs. Inx for a greater and less than the periodic value (1.9444) 
are quahtatively different. If a < 1.9444 the outer density slope increases without hmit. 
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while ii a > 1.9444 there is a maximum value for 7, followed by oscillations. 

If we abandon the 7 i„ ^ initial condition and instead use 7 ^„ > we can "overpower" 
the positive damping term and make 7 attain large values. The solution obtained with 
Jin — 1-05 and a — 1.975 is plotted as the long-dash line in Figure 14. The price we pay for 
resorting to large 7^„ is non-zero dln{pa'^)/dlar at the center of the halo. (Furthermore, 7 
decreases towards the center, and becomes negative: these halos have density holes in the 
center.) It is impossible to make a halo with a > 1.9444, constant pressure core, and outer 
halo 7 > 3 — 4 within the confines of the spherically symmetric, isotropic velocity dispersion 
Jeans equation. 

3.5.1. NFW vs. Hernquist-type density profiles 

We illustrate the positive and negative damping solutions with Hernquist and NFW 
density profiles, respectively. The recent work of Kazantzidis ct al. (2004) provides us with 
N-body realizations of Hernquist (their Model D) and NFW (Model F) profiles. These 
halos are not the result of the usual cosmological evolution starting from small density 
perturbations; instead they were explicitly constructed as isolated equilibrium objects. To 
qualify for the analysis using equation 6 the halos must be in equilibrium, have isotropic 
velocity dispersion ellipsoids, and power-law phase-space density distributions. Models D 
and F meet the first two requirements almost exactly, while the last requirement is met 
approximately, as we describe below. The N-body generated halos have been evolved for 
50 crossing times and thus have been explicitly shown to be in equilibrium. The velocity 
ellipsoid parameter, /3, equation 1, is zero for Models D and F. To estimate the phase-space 
density distribution we used Figure 3 of Kazantzidis ct al. (2004) to measure density and 
velocity dispersion profiles at seven points over 2-2.5 decades in radius. 

The phase-space density distribution of Model D is reasonably approximated by a power- 
law, but it is not an exact power-law: the slope based on the inner half of the halo is a = 2.01, 
while the outer half has a = 2.53. The average over the whole halo is with a = 2.27 ±0.11; 
the error is the RMS dispersion about the best fit line. We use this a value to generate 
the Jeans equation solution, and plot it as the dot-dash line in Figure 14. The solution 
and the actual data points from Kazantzidis et al. (2004) (empty squares) are in very good 
agreement. For reference, we also plot the Hernquist density profile, p oc [x{l + x)^]~^, as a 
thin solid line. 

The phase-space density profile of Model F is well fit by a single power-law with a = 
1.95 ± 0.04. This value is in excellent agreement with the periodic solution a — 1.9444 that 
we derived from the Jeans equation in Section 3.4. Since Kazantzidis et al. (2004) did not 
set out to produce halos with a close to the periodic value, we speculate that 1.9444 is an 
attractor, albeit a weak one: halos whose properties are already close to those of the periodic 
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solution (because of input constraints, dynamical evolution, or some other reason), end up 
with a 1.9444. The range of a values exhibited by the ESIM halos in Section 2 also attests 
to the weakness of the "periodic solution" attractor. 

The density profile of Model F (solid dots in Figure 14) can be reproduced with the 
Jeans equation solution with a — 1.95 and large 7^„, around 1, or with a smaller a and 

= 0. For example, a — 1.875 and 7i„ = solution is plotted as a thick sohd line; the 
NFW empirical fit (thin sohd line) is plotted for reference. 

While the Model F value of a = 1.95 ± 0.04 is fully consistent with the periodic solu- 
tion, it is inconsistent with TNOl's 1.875 at ~ 2(7. There are at least two possibilities why 
the a value obtained from N-body simulations (1.875) deviates from the periodic solution 
(1.9444) obtained assuming symmetric, isotropic halos in perfect hydrostatic equilibrium. 
N-body halos are subject to dynamical evolution, like gravitational collapse and mergers. 
Any of these could drive a away from the periodic solution value. (We remind the reader 
that Kazantzidis et al. models did not undergo the type of dynamical evolution that cos- 
mological halos undergo.) Alternatively, it is possible that our assumption of zero velocity 
dispersion anisotropy is not a good enough approximation. The actual anisotropy in halos 
has a gradient, with the inner region being isotropic and outer regions being more radial 
(Barnes et al. 2005a; Hansen & Moore 2004; Cole & Lacey 1996). Extending equation 6 
from /3 = to constant positive anisotropy /9 7^ 0, shifts the periodic solution to smaller 
a values. Specifically, a = 1.875 is the periodic solution for f3 ~ 0.3125 equilibrium halo. 
We postpone a detailed treatment of velocity anisotropy to a later paper. Finally, there are 
other, more subtle issues that can drive a away from its periodic value; halo triaxiality is 
one example. 

In this subsection we discussed Kazantzidis et al. (2004) models as examples of equilib- 
rium halos and compared them to the results of the Jeans equation analysis. However, these 
halos may not be of much help when trying to understand halos formed in cosmological sim- 
ulations, because of the very different dynamical evolution that the latter undergo. In other 
words, not all halos that are in equilibrium could have been formed in cosmological settings. 
A good example is provided by the Hernquist profile (Model D), which is in equilibrium, but 
does not describe halos generated in cosmological simulations. It follows that the values of 
a measured from the Hernquist profile are not particularly useful in elucidating the nature 
of halos evolved in cosmological settings, either N-body or ESIM. 



3.6. Conclusions from the Jeans equation analysis 

Our goals in Section 3 were two-fold: gain insight into the properties of solutions of the 
Jeans equation, and find out if the Jeans equation singles out any value of a. Analyzing 
equation 6 we found that the isotropic Jeans equation does indeed point to a special a: 
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a — 1.9444 results in periodic variation of 7 vs. \nx. Thus, a — 1.9444 divides the solutions 
into those where the damping term is typically negative, and 7 increases without limit, 
producing halos with density profiles that steepen at large radii, and those where the damping 
term is typically positive, and 7 asymptotes to a constant value, and hence power-law density 
profiles (provided the initial 7' is small). Interestingly, for the unique periodic solution 
a = 1.9444 the shallowest inner density slope and the steepest outer slope are 0.78 and 3.44 
respectively, values which are similar to those of N-body generated halos. 

For halos with constant non-zero velocity anisotropy, f3, the periodic solution shifts to 
smaller values of a, such that halos with purely radial orbits have periodic a = 1.7222. The 
N-body simulated halos have average anisotropics between and 1; so the corresponding 
periodic solution a should be bracketed by 1.7222 and 1.9444. Given this information, it is 
not surprising that the N-body generated halos that have undergone considerable dynamical 
evolution and have /3 > are characterized by an intermediate a value, a ~ 1.88. Our own 
tests with ESIM halos (Section 2) produce a range of a values, depending on the parameters 
of the initial set-up, and the properties of final halos. In a sense, ESIM halo results agree 
with the conclusions from the Jeans equation analysis because neither predict that all halos 
should have the same a. Based on the analysis of Section 2 and 3 we conclude that the radial 
slope of the phase-space density, a, is not fixed by the condition of hydrostatic equilibrium, 
and so must be the result of some other physics. 



4. Conclusions 

Examination of a handful of N-body dark matter halos indicates that the phase-space 
density, p/o"^, is a power-law in radius, even though the density profile is not (Taylor & 
Navarro 2001). To determine if the scale-free nature of p/(J'^ is peculiar to N-body halos 
(may be even just those presented in TNOl), or is a more generic outcome of collisionless 
gravitational collapse, we turned to a different method of generating equilibrium halos. The 
Extended Secondary Infall Model (ESIM) halo formation, which we used in Section 2, pro- 
ceeds very differently from that of N-body simulations. ESIM halos are spherically symmet- 
ric, and non-rotating. By assumption, the collapse conserves the angular and radial momenta 
of all shells individually. The final equilibrium halos are supported against collapse by radial 
and tangential velocity dispersion. 

Despite the differences between N-body and ESIM methods, and the resulting density 
profiles, both types of halos are characterized by a scale-free p/a^ profile. We speculate that 
because this feature is common to the two types of halos, it is an outcome of the one piece of 
physics shared by the two methods: violent relaxation. At the very least violent relaxation 
is a sufficient, if not the necessary condition needed to produce p/a^a: r~". We will further 
test this hypothesis in subsequent work. 
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If the scale-free nature of p/cr^ is, in fact, a generic feature of coUisionless halos, it implies 
an additional constraint on the halo structure as the density and velocity dispersion profiles 
are now closely related. In particular, the analysis of the Jeans equation and subsequent 
implications for halo formation are now not as unconstrained as before. Using this additional 
condition, we have rewritten the second derivative of the Jeans equation as equation 6, 
and used it as the starting point of our analysis in Section 3. This analysis yielded some 
interesting and novel insights into the solutions of the Jeans equation. In particular, we 
identified the unique periodic solution of equation 6, where the halo density profile slope 
oscillates between a minimum values of 0.78 and a maximum value of 3.44. These density 
profiles slopes are close to the values found in N-body simulated halos. The corresponding 
a — 1.9444 is the only special a value admitted by the spherically-symmetric, isotropic Jeans 
equation. 

In the future, we plan to extend the analysis of the Jeans equation to include non-zero 
velocity dispersion anisotropy distributions. Our ultimate goal is to identify the physical 
processes that shape the dark matter halos produced in simulations, and eventually, in the 
real Universe. 
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A. Convex density profile solutions of the Jeans equation 

All analytical density profiles that appear in the literature and describe mass (or light) 
distribution in galaxies are "concave", i.e. the double logarithmic density slope steepens 
with radius, 7' > 0. However, the Jeans equation also admits solutions with 7' < 0. These 
can be obtained by integrating equation 6 starting with 7j„ given by |[q; — 7] -|- 1 = and 
7^„ < 0. The density profiles asymptote to constant 7j„ = ck -1-1.5 at the center, and constant 
lout — (2/5)0; in the outer parts. In the representation of Figure 12 these profiles start at 
the upper thick solid line and end at the lower thick solid line. For any value of a the exact 
analytical solution is, 

I'^-l'm + lil- 70)', 70 = (14a + 15)/20, ^'^ = 3(2a + 5)7200, (Al) 
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and the double logarithmic density profile slope is given by 



7 = -{^iLm''' tanh [{2^'J^flHnx\-rio. (A2) 

The form of these equations is similar to that of equations 8-9 that describe concave density 
profiles (Section 3.4). The value of the density profile slope where the density slope changes 
most rapidly, 70, is depicted by a thin solid line in Figure 12. An example of this type of 
density profile is the lower half (7' < 0) of the long dash line in Figure 13. In these halos the 
velocity dispersion, cr, after reaching its minimum value at the radius where 7 = 70, starts 
to increase in the outer parts of the halo. The pressure, pcr^, asymptotes to a constant in 
the outer parts because 7out = (2/5)q; is constant there (see Section 3.2). 



B. Phase-space density of halos and the secondeiry infall model 

In Section 2 we examined the phase-space density of halos formed via the Extended 
Secondary Infall Model, ESIM, which differs from the original. Secondary Infall Model, SIM 
(Gunn & Gott 1972; Gott 1975) in many respects. The most notable differences are the 
shape of the pre-collapse peak (extended in ESIM, point mass in SIM), and the inclusion 
of secondary perturbations, and hence random motions in the ESIM (SIM halo collapse is 
purely radial). With these differences in mind, it is interesting to compare a values obtained 
in Section 2 to the prediction of SIM. 

In a spherically symmetric SIM halo collapse proceeds as follows. Initially, a small con- 
stant central mass excess, 5Mq is surrounded by material of average density. The dynamics 
of the pre-turn-around period is described by the parametric equations, r oc (1 — cos^), and 
t oc {9 — sin^). The maximum radius for a shell of initial comoving radius n is = ri5~^, 
where S is the initial fractional overdensity inside r^. Enclosed mass, M(< r^) oc prf, so 
Vi oc M^/^. And, S — SMq/M oc M~^. Combining these scahngs we get oc M^/^, for 
each shell in the halo. After reaching turn-around, a shell collapses back somewhat; it is 
usually assumed that the collapse factor, i.e. the ratio of the turn-around radius, r^, to the 
final radius in the virialized halo, is the same for all shells. Further assuming that each shell 
spends most of the time at its apocenter, which is proportional to r^,, one derives the result- 
ing density distribution in the final virialized halo: p(r) oc dM/dr, or, p oc M^^ oc Vm^^. 
This is a well known result. 

Let us take this approach a little further. In the real Universe the collapse of material 
will not be purely radial; there will be some random motion of particles, and associated 
kinetic energy. If we assume that during the collapse the kinetic energy is derived from the 
gravitational potential energy of infalling material, then the two types of shell energy are 
proportional, and oc GM/rm- Kinetic energy of random motion gives an estimate of the 
velocity dispersion: a oc {M/rmY^'^ ^ f^^^. Combining p(r) from the previous paragraph 
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with cr(r) , we get the radial dependence of the phase-space density, p/ oc (r~^/^) / {r"^/^) oc 
j,-i5/8^ i.e. a — 1.875. If the kinetic energy is not simply proportional to the potential energy, 
the final value of a will differ accordingly. The SIM phase-space density shape is always a 
power-law, by definition. 

These are the predictions of the highly simplified SIM. Given that ESIM halos of Sec- 
tion 2 have more complexity, and a variety of initial conditions, it is not surprising that ESIM 
halos show a wide range of a values. However, the typical values are not far from the SIM 
prediction of 1.875. Finally, the fact that a = 1.875 derived here, has the same numerical 
value as the best-fit a quoted by Taylor & Navarro (2001) for their N-body generated halos, 
is probably a coincidence, because the hierarchical halo formation characteristic of N-body 
simulations is qualitatively different from the highly idealized SIM. 

REFERENCES 

Avila- Reese, V., Firmani, C. & Hernndez, X. 1998, ApJ, 505, 37 

Bardeen, J.M., Bond,J.R., Kaiser, N. & Szalay, A.S. 1986, ApJ, 304, 15 

Barnes, E.I., Wilhams, L.L.R., Babul, A. & Dalcanton, J.J. 2005a, ApJ, accepted (astro- 
ph/0508160) 

Barnes, E.I., Wilhams, L.L.R., et al. 2005b, in preparation 
Bertschinger, E. 1985, ApJS, 58, 39 

Binney, J. & Tremaine, S. 1987, Galactic Dynamics (Princeton: Princeton University Press) 

Cole, S. & Lacey, C. 1996, MNRAS, 281, 716 

Dehnen, W. 2005, preprint, astro-ph/0504246 

Dehnen, W. & McLaughlin, D. 2005, preprint, astro-ph/0506528 

Diemand, J., Zemp, M., Moore, B., Stadel, J. & CaroUo, M. preprint, astro-ph/0504215 

Gott, J.R. 1975, ApJ, 201, 296 

Gunn, J.E. & Gott, J.R. 1972, ApJ, 176, 1 

Hansen, S. H. 2004, MNRAS, 352, L41 

Hansen, S.H., Egh, D., HoUenstein, L. & Salzmann, C. 2005, New Astron., 10, 379 
Hansen, S. H. & Moore, B. 2004, astro-ph/0411473 



22 



Huss, A., Jain, B. Steinmetz, M. 1999, ApJ, 517, 64 

Kazantzidis, S., Magorrian, J., Moore, B. 2004, ApJ, 601, 37 

Merritt, D., Navarro, J.F., Ludlow, A. & Jenkins, A. preprint, astro- ph/05025 15 

Moore, B., Governato F., Quinn, T., Stadel, J. & Lake, G. 1998, ApJ, 499, L5 

Moore, B., Quinn, T., Governato P., Stadel, J. & Lake, G. 1999, MNRAS, 310, 1147 

Navarro, J.F., Frenk, C.S. & White, S.D.M. 1997, ApJ, 490, 493 

Navarro, J.F. et al. 2004, MNRAS, 349, 1039 

Ryden, B. S., Gunn, J. E. 1987, ApJ, 318, 15 

Shapiro, P.R., Iliev, I.T., Martel, H., Ahn, K. & Alvarez, M.A. In Progress in Dark Matter 
Research, eds. Blain, J.V., (Nova Science Publishers, Inc, 2004) 

Syer, D. & White S. D. M. 1998, MNRAS, 293, 337 

Taylor, J.E. & Navarro, J.F. 2001, ApJ, 563, 483 

Tremaine, S., Richstone, D. O., Byun, Y.-L, Dressier, A., Faber, S. M.,Grillmair, C., Kor- 
mendy, J. & Lauer, T. R. 1994, AJ, 107, 634 

WiUiams, L.L.R., Austin, C., Barnes, E., Babul, A. & Dalcanton, J.J. Proceedings of 
" Baryons in Dark Matter Halos" . Novigrad, Croatia, 5-9 Oct 2004, eds: R. Dettmar, 
U. Klein, P. Salucci. Published by SISSA, Proceedings of Science, http://pos.sissa.it, 
p.20.1 

WiUiams, L.L.R., Babul, A. & Dalcanton, J.J. 2004, ApJ, 604, 18 



This preprint was prepared with the A AS macros v5.2. 



23 




-4 -2 Z -4 -3 a 

Fig. 1. — Density vs. radius in units of virial radius, r2oo, for galaxy-mass halos generated 
with Th = O.OSMpc, and perturbation amplitude scaling parameters of 0.4 [A], 0.6 [B], 1 [C], 
and 2 [D]. Density is expressed in units of the critical density, and multiplied by (r/r2oo)^^'S 
to accentuate deviations from a power-law. —Jfu, the best power-law fit to the density 
profile between 10~^'^ and 1 virial radii, is noted in each panel. The vertical normalization 
in this plot is arbitrary. The dotted, short dash and long dash hues represent radii where 
the density slope is 1, 2 and 3; 2,3. The NFW profile is over-plotted as a long dash line 
in all 4 panels. The NFW concentration parameters were computed using the masses of the 
corresponding halos. 
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Fig. 2. — The final properties of galaxy-mass halos with = O.OSMpc as as a function 
of the perturbation amplitude scahng parameter. [A]: Slope of the log-log density profile, 
jfit] [B]: Slope of the log-log p/a^ profile, afu] [C]: r^=i, r^=2, and r^=3: radn, in terms of 
virial, where the double logarithmic density slopes attain values of 1, 2, and 3, respectively; 
[D]: Virial mass of the halo (squares, and left vertical axis) and percentage mass lost during 
formation, with respect to the final halo mass (crosses, and right vertical axis). Note that 
as the amplitude of secondary perturbations is increased, the resulting halos have shallower 
density slopes ([A] and [C]) due to increased angular momentum, steeper phase-space density 
slopes [B] , and more mass loss [D] due to the increased probability of matter reaching escape 
velocity. 
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Fig. 3. — The top half of each of the four panels shows the RMS velocity profile for halos with 
perturbation amplitude scaling parameters of 0.4 [A], 0.6 [B], 1 [C], and 2 [D]. The lower half 
of each panel shows the velocity dispersion anisotropy, f3. Again, the values of r^=i,2,3 are 
shown in dotted, short dash, and long dash lines. Halos [A] and [B] have somewhat radial 
velocity anisotropics, while halos [C] and [D] have nearly isotropic velocity dispersions up to 
the virial radius. 
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Fig. 4. — The anisotropy parameter (3 vs. the density slope ^jfn relation for ESIM halos 
with scaling parameters of 0.4, 0.6, 1, and 2. Each set of points represents a halo; each point 
represents a different radial section of that halo. We also display the results of the 'group' 
and 'cluster' sized halos discussed in Section 2.2.5. 
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Fig. 5. — The phase-space density profiles of halos with perturbation amphtude scahng 
parameters of 0.4 [A], 0.6 [B], 1 [C], and 2 [D]. The p/a^ data have been normahzed to 
a power-law fit between 10~^'^ and 1 virial radius. The value of the power-law index fit is 
indicated in each panel as afu- Also, the slope of the p/a^ profiles from rjr^m of [10^^"'', 10^^] 
and [lO^*^'^, 1] arc noted as aiow and cihigh, for comparison. Note the nearly scale-free behavior 
for 2-3 decades in radius. 
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Fig. 6. — The hydrostatic equihbrium (HSE) condition: the LHS and RHS of equation 2 
are shown as points and the smooth curve, respectively. Panels A though D represent halos 
with scaling parameters of 0.4, 0.6, 1, and 2, respectively. Where the line traces the plotted 
points, the halo is in hydrostatic equilibrium. The dotted, short dash and long dash lines 
vertical lines are the same as in Figure 1. 
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Fig. 7. — Same as Figure 2, but for group-mass halos with rh — 0.35Mpc. 
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Fig. 8. — Same as Figure 2, but for cluster-mass halos with = O.TMpc. 
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log r/r200 log r/r200 

Fig. 9. — Profiles of two halos, both having scahng parameter of 1, generated with Vh — 
0.35Mpc (panels A through E) and Vh = 0.7Mpc (panels F through J). Shown, from top 
to bottom, are the density profiles, with the power- law best fit taken out, as in Figure 1 
(the NFW formula is over-plotted as a long dash line), the RMS velocity dispersion, the 
velocity anisotropy f3, the phase-space density with the power-law best fit taken out, and 
the equilibrium condition, showing the LHS (points) and the RHS (smooth curve) of eq. 2. 
Again, the values of jfit, afu, aiow, and ahigh are noted in the appropriate panels. Note that 
the phase-space density profiles are well fit by power-laws, while the density profiles are not. 
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Fig. 10. — The double logarithmic slope, afu, of p/a^ vs. radius, as a function of the virial 
halo mass, for three different P{k) filtering scales, = 0.05,0.35 and 0.7Mpc (triangles, 
squares, and pentagons). Filled (empty) symbols show afu values computed using 10~^'^ —>■ 
1 (lO^^'^ 10~°'^) radial range, given in terms of r2oo- Note that aju is not uniquely 
determined by halo mass, amplitude of secondary perturbations, halo filtering scale, etc. 
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Fig. 11. — An example of a halo with discontinuous perturbation amplitude scaling param- 
eters: shells with the initial comoving radii between 0.59 and 0.79 Mpc were assigned a 
scaling parameter of 5, while the rest of the halo (up to the last shell at 1.92 Mpc) were as- 
signed a scaling parameter of 0.9. [A]: Density profile, normalized to a power-law fit between 
of 10^^'^ and 1; [B]: RMS velocity dispersion profile in km/s; [C]: Velocity dispersion 
anisotropy profile, (3. [D]: Phase-space density profile, normalized to its power-law fit; [E]: 
The equilibrium condition of equation 2: the dots represent the LHS and the thin line the 
RHS. Where the hne traces the plotted points, the halo is in hydrostatic equilibrium. 
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Fig. 12. — Pictorial depiction of equation 6. Thick solid and long-dash hnes represent the 
terms on the LHS of equation 6. The dot-dash line is the coefficient of the damping term 
on the RHS of equation 6. The two thin solid lines are 70(a) 's given in equations 8 and 
Al. The three vertical dotted lines at a = 1.875, 1.9444 and 2 arc plotted to guide the eye. 
Periodic solutions arc obtained for a = 1.9444 (Section 3.4). A set of closely spaced crosses 
are afu and jfit values for galaxy-mass ESIM halos with /3 ~ 0. 
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Fig. 13. — 7' vs. 7 for a few Jeans equation solutions. (7 is the log-log halo density 
profile slope.) The solid line approximates the NFW fit. The long-dash line is the limit 
cycle for equation 6: a = 1.9444, initial 7j„ = (2/5)q; + 0.001 = 0.7788, initial 7'^^ = 0. 
The dotted line (closed loop) also has a = 1.9444, but different initial = 2.7111, and 
7-^ = 0. The solid {a = 1.875) and short-dash {a — 1.975) lines are the same as in Figure 14. 
For both, 7i„ = (2/5)q; and 7 •„>0. The top arc (parabola) of the long-dash line has an 
analytic description, given by equations 8 and 9; the bottom arc (parabola)is described by 
equations Al and A2. 
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Fig. 14. — Double logarithmic density profile slope, 7 vs. radius, x (arbitrary normalization). 
The thick solid line is the solution of the Jeans equation that best approximates the NFW fit: 

a. = 1.875, 7 •„>0. Short-dash (a = 1.975, 7 •„>0), long-dash {a = 1.975, 7'^ = 1.05) and 
dot-dash (a = 2.27, 7'^^^ = 2.56) lines are other Jeans equation solutions. All four solutions 
have Jin = (2/5)0;; sec Sections 3.4 and 3.5 for details. The solid and short-dash lines are 
the same as in Figure 13. The empty squares and solid dots represent Models D (Hcrnquist) 
and F (NFW) of Kazantzidis et al. (2004). The two thin solid lines are the Hernquist (left) 
and NFW (right) density profile fits. 
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